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Abstract 

We present evidence from computer simulations for glassy dynamics in suspensions of monodis- 
perse hard ellipsoids. In equilibrium, almost spherical ellipsoids show a first order transition from 
an isotropic phase to a rotator phase. When overcompressing the isotropic phase into the rota- 
tor regime, we observe super- Arrhenius slowing down of diffusion and relaxation, accompanied by 
two-step relaxation in positional and orientational correlators. The effects are strong enough for 
asymptotic laws of mode-coupling theory to apply. Glassy dynamics are unusual in monodisperse 
systems. Typically, polydispersity in size or a mixture of particle species is prerequisite to pre- 
vent crystallization. Here, we show that a slight particle anisometry acts as a sufficient source of 
disorder. This sheds new light on the question of which ingredients are required for glass formation. 

PACS numbers: 64.70.pv, 61.20.Ja, 61.25.Em, 82.20.Wt 
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Hard-particle models play a key role in statistical mechanics. They are conceptually 
and computationally simple, and they offer insight into systems in which particle shape is 
important, including atomic, molecular, colloidal, and granular systems. Ellipsoids are a 
classic model of non-spherical particles. We report here that this simple anisometry can 
hinder crystallization and facilitate glassy dynamics - a phenomenon which does usually 
not occur in monodisperse systems. Typically, polydispersity, additional particle species, or 
other sources of disorder are needed for the development of glass-like behavior, i.e. drastic 
slowing down of the dynamics without a change in structure. 

In recent years, there have been two studies closely related to our topic, which we briefly 
summarize here: Letz and coworkers [11] have applied idealized molecular mode-coupling 
theory (MMCT H, [sl) to the hard-ellipsoid fluid. In addition to conventional mode-coupling 
theory (MCT) j4|, MMCT takes orientational degrees of freedom into account. For nearly 
spherical ellipsoids, they predicted a discontinuous glass transition in positional and ori- 
entational degrees of freedom. MMCT locates the transition inside the coexistence region 
between the isotropic fluid and the positionally ordered phases. In addition, a continuous 
transition was predicted upon further compression into the rotator regime. This transition 
affects only the odd-parity orientational correlators, e.g. 180° flips. 

However, MCT cannot make statements about crystal nucleation. Hence, the MCT 
prediction of a glass transition is not sufficient to conclude that the transition will occur in 
a simulation or experiment. A prominent example is the overcompression of monodisperse 
hard spheres. Here, the nucleation barrier can be easily crossed, and crystallization always 
prevents glass formation. 

De Michele et al. jsl have recently studied the dynamics of hard ellipsoids by molec- 
ular dynamics simulations. The states which they simulated were mostly located in the 
isotropic region. They computed isodiffusivity lines, which showed that the dynamics of the 
positional and orientational degrees of freedom were decoupled, since the positional isodif- 
fusivity lines crossed the orientational ones at nearly 90°. This decoupling also appeared in 
correlation functions. The self-part of the intermediate scattering function displayed slight 
stretching only when overcompressing nearly spherical ellipsoids, while the second-order ori- 
entational correlator showed such stretching only for sufficiently elongated particles, i.e. near 
the isotropic-nematic transition. Clear indicators of glassy dynamics, however, would include 
a strong increase of relaxation times with volume fraction, even pointing towards dynamical 
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arrest. Typically, correlators then develop a two-step decay, whose second step is affected 
by this slowing down. Such phenomena were not seen in |5| as overcompression was not 
significant. 

We have performed Monte Carlo (MC) and molecular dynamics (MD) simulations for 
hard symmetrical ellipsoids of length-to- width ratios l/w = 1.25 (prolate) and l/w = 0.8 
(oblate). We overcompressed these systems into the rotator regime (i.e. beyond volume 
fraction (f) 0.55 6j). We found two-step relaxation both in positional and even-parity 
orientational correlators. Positional and orientational relaxation slow down more strongly 
than an Arrhenius law. Odd-parity orientational correlators indicate that flipping is not 
affected. The observed glassy dynamics are strong enough to compare with MMCT. Also, 
we compare the MC and MD results. 

The systems were equilibrated using MC at constant particle number A^, pressure P and 
temperature T 7|, |8| in a cubic box with periodic boundaries. Each system contained more 
than 3000 particles. Random isotropic configurations were used as starting configurations. 
Towards the end of each run, a configuration with a volume fraction close to the average 
volume fraction was chosen and scaled to the average volume fraction exactly. The systems 
were considered equilibrated 



23l | when the volume fraction had settled and all positional and 
the orientational correlators were independent of absolute simulation time. For production, 
we used MC and MD. In the MC simulations (now at constant volume) the step sizes were 
fixed to small values, so that unphysical grazing moves were negligible. The particles then 
mimic Brownian motion, similar to colloidal ellipsoids suspended in a liquid. The step 
sizes were the same for all runs, hence a unique time scale could be established. The MD 

CI 

simulations implemented free flight and elastic collisions [9|. 

The systems had a slight tendency to crystallize to the rotator phase at high volume 
fractions. To monitor crystallinity, we computed the local positional order using the bond- 
orientational order parameter lOj. The fraction of particles which were part of crystalline 
clusters never exceeded 2.6% and was typically below 0.5%. 

To demonstrate the slowing down of the positional degrees of freedom, we consider the self 
part of the intermediate scattering function, Fs{q,t) = (exp[zq- Ar(t)]), where q is the wave 
vector, Ar(t) the displacement of an ellipsoid after time t, and the angle brackets denote 
average over particles and ensemble average. In isotropic systems, is a function of the 
absolute value q only. Decay of Fs{q,t) indicates that structural relaxation has occured on 
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FIG. 1: Self-intermediate scattering functions at several volume fractions (f) (prolate: 0.474, 0.511, 
0.533, 0.551, 0.565, 0.578, 0.588, 0.598; oblate: 0.473, 0.504, 0.533, 0.550, 0.565, 0.578, 0.589, 0.598, 
0.606, 0.614). At high volume fractions there is a plateau on intermediate time scales. The final 
relaxation is slowed down strongly with increasing volume fraction, indicative of glassy dynamics. 



the length scale set by q. In Fig. [T] we show Fs{q, t) for all simulations. The wave length qma,x 
was chosen close to the first maximum of the static structure factor (prolate: 6.85/ w, oblate: 
7.85/w), i.e. it corresponds to the neighbor spacing. One can clearly see the development of 
a plateau with increasing volume fraction. This means that there are two distinct stages of 
relaxation, and the latter slows down dramatically upon increase of volume fraction (note 
ihe logarithmic time scale). This two-step decay is a typical phenomenon in glass formers 
^. It is interpreted in terms of particles being trapped in cages formed by their nearest 



4 




10^ 10^ 10^ 10^^ 10^ 10^ 10^ 10^ 10^^ 10^ 
MC steps MD steps 



FIG. 2: (Color online) Second-order orientational correlators at several volume fractions (p as in 
Fig. [H Again a plateau develops with increasing (j). Hence, orientational degrees of freedom are 
coupled to the positional ones. Also shown are the third-order correlators at the highest density 
(bold green lines). They do not slow down, indicating that flipping modes are not affected. 



neighbors. The initial decay corresponds to motion within the traps, and the final decay to 
escape. At high volume fractions, exceedingly cooperative rearrangements are required for 
escape, making such events rare. Indicators of caging have not been seen in hard ellipsoids 
before (unless the moment of inertia was strongly increased 111]). 

Dynamics in the orientational degrees of freedom is observed in terms of the second-order 
orientational correlator L2{t) = (l/2)(3 cos^ 0(t) — 1), where 6{t) is the angle between the 
orientation at time t and the original orientation of an ellipsoid. Decay of this function 
indicates that relaxation of orientational degrees of freedom has taken place. Since L2(t) is 
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an even function in cos 0{t), the head-to-tail symmetry of the particles is taken into account. 
Figure [2] shows the orientational correlation functions. As in the intermediate scattering 
functions, plateaus develop at high volume fractions. Evidently, the cages hinder rotations 
of the ellipsoids. As a consequence, orientational and positional degrees of freedom are 
coupled. This is in contrast with the decoupling found at lower volume fractions [5]. 

The shape of both positional and orientational correlators differs between MC and MD 
on short time scales, reflecting the individual microscopic dynamics. On intermediate and 
long time scales, they do not differ significantly. Furthermore, when the correlators of the 
highest few volume fractions are rescaled by their decay time, their long time parts fall onto 
a master curve. These properties confirm predictions of MOT 

Unlike L2{t), the third-order orientational correlation function L^lt) = (1/2) (5 cos^ 9{t) — 
3cos6{t)) does not show plateaus (bold green lines in Fig. [2]). Hence, while the overall 
reorientation slows down, flipping is barely hindered. This is in accord with the MMCT 
prediction of Letz et al. j| and has also been found for the case of diatomic Lennard- Jones 
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dumbbells [13[, and symmetric Lennard- Jones dumbbells [iJ]. We note in passing that 
crystallization, if it occurs, releases the orientational degrees of freedom: The orientational 
correlators accelerate by three orders of magnitude and no longer have a plateau. 

Next, we show that the slowing down in our systems is strong enough to test MOT 
asymptotic laws. To this end, we first return to the intermediate scattering functions and 
focus on their g-dependence at high volume fractions. Figure [3] presents these functions for 
a) the oblate system (MC) and b) the prolate system (MD; in each case the other method 
shows similar results), with wave lengths in the range 2.8 < qw < 20. We test two MCT 
predictions for the vicinity of the glass transition ^. Firstly, for the late stages of the 
plateau and the early stages of the final decay, these functions should be well-fitted by 
the von-Schweidler law (incl. the second-order correction), Fs{q,t) = /q — hqh^ + h'^h'^'', 
where /q is the plateau height, hq^ are amplitudes, and 6 is a system-universal exponent 
(also independent of the microscopic dynamics). Agreement is excellent, as shown for two 
examples in Fig. [3] (dashed lines). MC and MD data of both systems are consistent with 
b = 0.65 ± 0.2. Secondly, where the Fs{q,t) are near their plateaus, they should obey 
^s{(l,t) = fq + hqG(t), where hq is an amplitude, and G(t) is a system-universal function. 
This relation entails the "factorization property", i.e. that Fs{q,t) can be factorized into a 
g-dependent and a t-dependent part. To test this property for our systems, we transform 
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FIG. 3: Self-intermediate scattering functions Fs{q,t) of a) the oblate system (MC data) and b) 
the prolate system (MD data) at the highest volume fraction, for the range 2.8 < qw < 20. The 
dashed lines show fits to the von- Schweidler law Fs{q,t) = /q - h^qh^ + h^^h^^ with b = 0.65. c) 
and d) show the same data transformed to R{t) = [Fs{q,t) — Fs{q,ti)]/[Fs{q,t2) — Fs{q,ti)]. The 
collapse of the functions onto master curves demonstrates the factorization property. 



Fs{q,t) to R{t) = [Fs{q,t) - Fs{q,ti)]/[Fs{q,t2) - Fs{q,ti)], as done in [l5|, where h and 
t2 are times in the regime where the property holds. Since R{t) is not a function of q, 
all correlators should fall onto a single master curve. Moreover, the curves should remain 
ordered, i.e. a curve which is above another on the left-hand side remains above the other on 
the right-hand side. Panels c) and d) of Fig. [3] demonstrate the validity of the factorization 
property, and indeed they remain ordered. 

We finally demonstrate the applicability of an MCT scaling law for the slowing down, 
which even allows us to extract the MCT glass transition volume fraction. Fig. H] shows 
relaxation times and diffusion constants as a function of volume fraction. The relaxation 
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FIG. 4: Upper panels: Inverse relaxation times obtained from -Fs(9max,i) (triangles) and from 
L2{t) (circles), and diffusion constants D (squares), as a function of volume fraction MC data 
(filled symbols, lines to guide the eye, and rescaled to match MD time scale) and MD data (open 
symbols) show excellent agreement. Lower panels (MC data only): t~^^'^ multiplied by arbitrary 
factors for clarity, from Fs{q,t) (several values of q) and L2{t), demonstrating the validity of the 
MCT scaling law oc (^c — (t^V- 7 = 2.3. The straight-line fits indicate a glass transition volume 
fraction of (f)c = 0.615 it 0.005 (prolate system) and (f)c = 0.618 ± 0.005 (oblate system). MD data 
agree. 
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times are defined as the times at which the positional correlators have reached the value 
0.1, and the orientational correlators have reached the value 0.02 (since the plateau is quite 
low). The diffusion constants were determined from the mean squared displacements via the 
Einstein relation D = lim^^oo J^:((^r(t))^)/6. The upper panels of Fig. Hlshow inverse relax- 
ation times r^^ obtained from Fs{qma,x,t) (triangles) and from L2{t) (circles), and diffusion 
constants D (squares). The MC values (filled symbols) have been multiplied by a factor 
of 24 (prolate) and 18 (oblate) to match the MD time scale (L2: 19 and 16, respectively). 
The errors are of symbol size or less. Note that in each system, one common factor for 
the positional variables yields excellent agreement between MC and MD (open symbols), 
in agreement with MCT. The factor for the orientational relaxation times need not be the 
same since it depends on the choice of the orientational MC move size. The slowing down of 
all variables is super- Arrhenius. According to MCT, it should approach a power law of the 
form D oc t^^ oc (0^ — 0)'^, where 0^ is the MCT glass-transition volume fraction, and 7 is 
related to the von-Schweidler exponent b. Both 0c and 7 should be system-universal. In the 
lower panels of Fig. H] we demonstrate the validity of this prediction for the MC results of 
Fs{q,t) (from top: prolate: qw = 6.85, 11, 16; oblate: qw = 7.85, 16, 20) and L2 (bottom). 
The exponent 7 = 2.3 was chosen in agreement with b = 0.65. The straight-line fits comply 
with a common 0c = 0.615 ± 0.005 for the prolate system and with 0c = 0.618 ± 0.005 for 
the oblate system. The fact that there is a common value for positional and orientational 
relaxation times further demonstrates the strong coupling of these degrees of freedom. We 
found agreement with the analogous analysis of the MD data. The values we found differ 
from the MMCT predictions of Letz et al. [1], viz. 0c = 0.540 and 0.536 {l/w = 1.25 and 
0.80, respectively). The mismatch between the numerical MCT calculations based on static 
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structure, and scaling law fits based on simulated dynamics is, however, not unusual 
and has been attributed to activated ( "hopping" ) processes for which MCT does not account. 
A similar mismatch is found in the hard-sphere system ISj. We note that the present study 
also displays the prolate-oblate symmetry seen in previous work on the equilibrium proper- 
ties and dynamics of ellipsoids [l|, 0|. However, we observed that crystallization does not 
have this symmetry: the prolate system crystallizes more readily. 

In summary, we have performed molecular dynamics and Monte Carlo simulations of 
the hard-ellipsoid fluid. In this very simple anisometric model we observe glassy dynamics 
sufficiently strong that MCT asymptotic scaling laws can be tested and are found to apply. 
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We find strong coupling of positional and orientational degrees of freedom, leading to a 
common value for the glass-transition volume fraction 0c for positional and orientational 
relaxation times {l/w = 1.25: 0^ = 0.615 ± 0.005, l/w = 0.80: 0^ = 0.618 ± 0.005). The 
presence of glassy dynamics has been predicted by MMCT. However, as MMCT cannot make 
a statement about crystallization, a test by simulation was required. We argue that particle 
anisometry acts as a sufficient source of disorder to prevent crystallization. This sheds new 
light on the question of which ingredients are required for glass formation. Experimental 
studies of glassy dynamics in the isotropic phase of liquid crystals have been conducted 
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[19[, but not in ellipsoids. It is possible to synthesize ellipsoids of colloidal size 
and to study their dynamics with confocal microscopy 22j. In the light of the above, an 
experimental study of glassy dynamics in colloidal hard ellipsoids seems very promising. 
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